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Abstract 



We derive a spectral representation for the oblate spheroidal wave operator which 
is holomorphic in the aspherical parameter in a neighborhood of the real line. For 
real f2, estimates are derived for all eigenvalue gaps uniformly in O. 

The proof of the gap estimates is based on detailed estimates for complex solutions 
of the Riccati equation. The spectral representation for complex Q is obtained using 
the theory of slightly non-self adjoint perturbations. 

> 

o . 

1 Introduction 

ID 

Recently an integral representation was derived for solutions of the scalar wave equation in 
the Kerr black hole geometry [I]. This result relies crucially on a spectral representation 
for the oblate spheroidal wave operator for complex values of the aspherical parameter f2 
(also referred to as "ellipticity parameter" or "semifocal distance"). In the present paper, 
this spectral representation is proved. The reason why this problem deserves to be worked 
out in a separate paper is that most of our methods apply in a much more general context. 
Namely, the core of the paper is to derive estimates for the eigenvalue gaps A n+ i — \ n for 
real which are uniform in 0, and n. To this end, we need to control the eigenvalues 
and the behavior of the wave functions in detail. Our method is based on invariant region 
estimates for the complex Riccati equation and applies to general Sturm-Liouville or one- 
dimensional Schrodinger problems. In particular, it gives refined error estimates for WKB 



approximations. We regard the spheroidal wave equation as a model problem for working 
out these estimates. 

Despite the vast literature on spectral estimates for the Schrodinger equation (see 
e.g. 12 and the references therein), gap estimates are rarely found in the standard liter- 
ature. Most papers are concerned with the two lowest eigenvalues (HI El, or they apply 
in special situations like for a a nearly constant potential |1U| . Probably, this is because 
gap estimates depend sensitively on the detailed form of the potential (as one sees in the 
example of a double- well potential), making it difficult to get general results. Our method 
requires that the potential is piecewise monotone and that we have good control of its 
derivatives. 

We now introduce our problem and state our results. The spheroidal wave equation 
is the eigenvalue equation for the spheroidal wave operator, a linear elliptic operator 
with smooth coefficients on the unit sphere S 2 . Since the spheroidal wave operator is 
axisymmetric, we can choose angular variables •& € (0, ir) and ip £ [0, 2ir) (with i? the angle 
to the axis of symmetry) and separate out the 92-dependence with the plane wave ansatz 
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4>{'d,ip) = e ^ 0(i?), k G Z. After this separation, the spheroidal wave equation takes the 
form 

^9 = AG, (1.1) 
where A is the linear differential operator of second order 

A = --^ S m 2 ti-^ + -^-(nsm 2 V + k) 2 (1.2) 
a cos v a cos v sin v 

on the interval t? G (0, 7r). Here fi £ C is the aspherical parameter. In the special 
case O = 0, the spheroidal wave operator simplifies to the spherical Laplacian, and the 
Legendre polynomials -/-^(cosi?) are explicit solutions to 1)1. lj) . We shall consider the 
spheroidal wave equation for fixed k, but for a variable complex parameter f2. The fact 
that the eigenfunction 4> should be smooth at the poles i9 = 0, it of the sphere gives rise 
to the following boundary conditions, 

lim e'(#) = iffc = 

lim 6(0) = iffc^O. 1 ; 

We consider *4 as an operator in the Hilbert space TL = L 2 ((0, n), sin i? di?) with domain of 
definition given by those functions in C 2 (0, tt) which satisfy the boundary conditions (|1.3|) . 
Note that the potential in the spheroidal wave operator is in general complex, 

Im / (Osin 2 ^ + fc) 2 \ = 2 (R B n s in 2 d + k)lmn, (1.4) 
\ sin "9 J 

and therefore A is symmetric only if fi is real. In previous works, asymptotic expansions 
for individual eigenvalues are derived > an d it is shown numerically that eigenvalues 
can degenerate for non-real O [7], but rigorous estimates or completeness statements are 
not given. Our main result is the following spectral representation for O in a neighborhood 
of the real line. 

Theorem 1.1 For any k G Z and c > 0, we define the open set U C C by the condition 

1 1 l + |Refi| J 

Then there is a positive integer N and a family of operators Q n (Q) on Ti defined for 
n G N U {0} and O G [/ with the following properties: 

(i) The Q n are holomorphic in f2. 

(ii) Qo is a projector on an N -dimensional invariant subspace of A. For n > 0, the Q n 

are projectors on one- dimensional eigenspaces of A with corresponding eigenvalues 
A n (f2). These eigenvalues satisfy a bound of the form 

\X n (0)\ < C(n)(l + \n\) (1.6) 

for suitable constants C(n). Furthermore, there is a parameter e > such that for 
flllnGN and O G U, 

|A n (0)| > ne. (1.7) 
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(iii) The Q n are complete, i.e. 

oo 
fc=0 

with strong convergence of the series. 

(iv) The Q n are uniformly bounded, i.e. for ollnGNU {0}, 

||Q n || < ci (1.8) 

with c\ independent of Q and k. 

If c is sufficiently small, c < 5, or the real part o/fi is sufficiently large, |Refl| > C(c), one 
can choose N = 1, i.e. A has a purely discrete spectrum consisting of simple eigenvalues. 

To avoid misunderstandings, we point out that by a "projector on an invariant subspace 
of A" we mean an operator Q which is idempotent and commutes with A. But Q will in 
general not be symmetric. 

In our proof we shall treat the imaginary part of the potential (|1.4j) as a slightly non- 
selfadjoint perturbation in the spirit of [HI V.4.5], see also [21 Chapter 12]. For this method 
to be applicable, we need good control of the eigenvalues of the corresponding selfadjoint 
problem. Our starting point is the following spectral decomposition of A in the case of 
real £1. 

Theorem 1.2 For any k 6 Z and !] £ K, the operator A has a unique selfadjoint ex- 
tension compatible with the boundary conditions 11. This extension, which we again 
denote by A, is a positive operator with compact resolvent and simple eigenvalues. It is 
invariant on the even and odd parity subspaces Ti^ defined by 

n = n + @n- with = {0 e n with ^ -■&) = ±<p($)}. 

We denote the eigenvalues of A restricted to Ti^ 1 by and count them with multiplicities, 

< xf < xf < xf < ■■■■ 

Using abstract methods (see [SJ Theorem 3.9, VII.3.5]), one could show that each 
eigenvalue A^(f2) has a holomorphic continuation to a neighborhood of the real axis. 
However, as pointed out in jSJ Remark 3.9, VII.3.5], this neighborhood will depend on n, 
making it impossible construct a neighborhood in which all the X^(Q) exist. Therefore, 
abstract methods only seem to give results which are much weaker than Theorem 11.11 
where the whole spectral decomposition is shown to have a holomorphic continuation to 
a neighborhood of the real axis. Furthermore, we point out that the parameter c in the 
statement of Theorem 11.11 can be chosen arbitrarily large. Therefore, the holomorphic 
family of operators Q n (£l) is not only defined in a small neighborhood of the real axis, but 
in a strip (|1,5[) which can enclose any bounded subset of the complex plane. The key for 
getting this strong result are the following gap estimates uniform in n and 0. 

Theorem 1.3 For any k 6 Z and 7 > 0, there is a positive integer N such that 

A n+i - A n > 7 f° r alln>N and !leK. 

7/7 is sufficiently small or \Q\ is sufficiently large, one can choose N = 1. 
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The paper is organized as follows. In Section [2] we prove Theorem 11.21 and reduce 
Theorem ll.3l to gap estimates for a self-adjoint Sturm-Liouville operator on the interval $ 6 
(0, 7^] with suitable boundary conditions. In Sections EHf3 we introduce the complex Riccati 
equation and develop general techniques for analyzing its solutions. In Section [7| and |H] we 
apply these techniques to the spheroidal wave operator and prove Theorem 11.31 Finally, 
in Section El we use perturbative methods to proof Theorem ll.il 



2 Basic Considerations 

Until the end of Section we will consider the spheroidal wave equation Ql.lj) for real f2. 
Using that (|1.2|) is invariant under the transformations f2 — > — f2 and k — > —k, we can 
assume throughout that 

n > o. 

Let us derive a spectral representation of the spheroidal wave operator. One possible 
method would be to apply elliptic theory to the spheroidal wave operator on S 2 before sep- 
aration of variables. After choosing a self-adjoint extension on the Hilbert space L 2 (S 2 ), 
one could apply the abstract spectral theorem, and projecting the resulting smooth eigen- 
functions on the subspace for fixed k would give the desired spectral decomposition for 
the ordinary differential operator (|1.2j) . For clarity, we will in this paper restrict attention 
to ODE techniques. Thus we avoid elliptic theory and prefer to apply Sturm-Liouville 
theory. In the variable u = $ G (0,7r), the operator (|1.2j) can be written as 

. Id. d 1 , . n , n2 

A = — Sinn- — | T) — (S2sm u + k) . 

sin u du du sin u 

In order to bring this operator to the standard Sturm-Liouville form, we introduce the 
function Y by 

Y = Vsin-u© . (2.1) 



Then 
where 



BY = XY 



1 d . d 1 1 . ,, , : 

B = : — sin u — : H 5 — ( s I sin u + k) 

Vsinu du du Vsinu sin z u 



d 2 1 cos 2 n r . / 1 \ 1 . . 2 .r 

—r-z + - — n vsinu 7= H 5— (S 2 sin u + ky 

du A 2 sm z u V v sin u J sin u 

d 2 1 COS 2 U 1 1/^.9 \2 

~r~o ~ t ~. — 9 + — p (£2 sin u + k) . 

du 2 4 sin u 2 sin u 



Thus y satisfies the Sturm-Liouville equation 

d 2 



where V is the potential 



du2+ V)F = 0, (2.2) 



U = O 2 sin 2 u + f k 2 - - ) —73 u (2.3) 

5 4/ sin 2 n 
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with n the constant 

H = X-2Qk + ^. (2.4) 

The transformation ()2.1|) from G to Y becomes a unitary transformation if the integration 
measure in the corresponding Hilbert spaces is transformed from sin u du to du. Thus the 
eigenvalue problem on 7i is equivalent to ()2.2|) on the Hilbert space 
The boundary conditions 1)1. 3|) at u = (and similarly at u = ir) can be written as 

Hm(ir5 Y)\u) = if = 

«\° (2 5) 

limu"5 y( u ) = if fc / . 

The Sturm-Liouville equation Q2.2JI is singular at the two end points u = 0, it. An 
asymptotic expansion near u = shows that (|2.2|) has fundamental solutions Y\h of the 
following form, 

Y 1 = y /u + 0(u*) , Y 2 = y/u log(u) + 0{vh) if As = 

y x = u l+l fe l + o( u l+l fc l) , y 2 = -u5-l fc l + o(nl-l fc l) iffc^O, 

and similary at u = 7r. In the case k ^ 0, Y\ is square integrable near u = 0, whereas y 2 is 
not. Thus, using Weyl's notation, the Sturm-Liouville operator is in the limit point case at 
both end points, and thus A is essentially selfadjoint (see |H Sections 9.2, 9.3] or Chapter 
XIII. 2]). In the case k = 0, on the other hand, both fundamental solutions are square 
integrable. This is the limit circle case, and the von-Neumann boundary conditions (|2.5|) 
choose a unique self-adjoint extension (see [21 Sections 9.4] or Chapter XIII.2]). We 
conclude that the Sturm-Liouville operator in (|2.2|) has a unique self-adjoint extension 
in -L 2 ((0, 7r)) which satisfies the boundary conditions (|2.5j> . Hence the spectral theorem for 
unbounded operators in Hilbert spaces gives us the desired spectral representation of A. 

For each A € M, there are (up to a constant) unique solutions of the ODE which satisfy 
the boundary conditions at u = and u = vr, respectively. If the Wronskian of these 
two solutions vanishes, we obtain an eigenfunction in L 2 ((0,7r)). Otherwise, these two 
solutions can be used to define the resolvent (see [31 XIII.3]), which is compact (see 
XIII.4]). This shows that the operator A has a purely discrete spectrum consisting of 
simple eigenvalues without limit points. The positivity of A is obvious from ([1.2)) , 

Finally, the boundary value problem ()2.2l I2.5j) is invariant under the transforma- 
tion u — > 7r — u. Hence the parity subspaces := {(p G L 2 ((0,7r)) | (f>{ir — u) = ±<j)(u)} 
are invariant under A. This concludes the proof of Theorem ll.2l 

Clearly, the eigenfunctions Y^ of even and odd parity satisfy at u = ^ the boundary 
conditions 

(y„ + )'(f)=o 

(A (2 - 7) 

Y~ - =0. 



.2- 

This makes it possible to consider instead of the the interval (0, ir) only the interval (0, ^] 
together with the boundary conditions (|2.5U2.7)) . In what follows, we shall always consider 
the boundary value problem (J2.21 12.51 12.7)1 . 

In order to better understand Theorem 11.31 it is useful to consider the limits n — > oo 
and £1 — > cxd. For fixed f2 and large n, Weyl's asymptotics applies and yields that the 
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eigenvalues of A behave for large n like the eigenvalues of the operator —-jr?, 

An ~ n2 and A^ +1 -A* ~ n. 

Therefore, it is obvious that the statement of Theorem 11.31 holds for any fixed O and 
sufficiently large N = N(Q). The estimate 

|A n (o)-A n (o')| < \\A{n) -A(pf)\\oo < |n-n'| (n + n' + 2|fc|) (2.8) 

yields that eigenvalues of A are locally Lipschitz in f2, uniformly in n. This shows that 
the constant iV(O) can be chosen locally uniformly in $7. If conversely we fix n, the nth 
spheroidal eigenvalue A n has for large O the asymptotic expansion (see [H] or |llj ) 

f 2(n + l)Q + 0(1) if n - k is even, 
A„(fi) = { (2.9) 
{ 2nQ + O(l) if n - k is odd. 

Hence for each n, we can make the eigenvalue gap arbitrarily large by choosing f2 suffi- 
ciently large. We conclude that it remains to show that the eigenvalue gaps are bounded 
uniformly as both iV and |fi| become large. This is the hard part of Theorem II .31 and we 
state it as a separate Lemma. 

Lemma 2.1 For any given k £ Z and c > 0, there are constants N e N and H,q > such 
that 

A^ +1 — A^ > c for all n > N and £1 > Qq. 

The proof of this lemma requires detailed eigenvalue estimates. We will complete it in 
Section [7J and this will also finish the proof of Theorem II .31 

Finally, the node theorem |15| Theorem 14.10] tells us about the number of zeros of 
the spheroidal wave functions. In our setting, the statement of the node theorem can 
easily be derived as follows. Using the initial conditions (|2.5|) together with (|2.8|) . we 
obtain from the Picard-Lindelof theorem that the eigenfunctions corresponding to the 
eigenvalue A^ depend smoothly on the parameter Q. Furthermore, the asymptotics near 
the boundaries (|2.51 12 .7(1 shows that the functions Y^ 1 have no zeros in the intervals (0, e) 
and — e, §), for a parameter e > which depends continuously on S7. Thus, if the 
number of zeros of the function Y^ changed at some Q, there would be a u £ (0, 7r) with 
Y^{u) = = (Y^)'(u), in contradiction to the fact that Y^ does not vanish identically. 
We conclude that the number of zeros of Y^ on (0, ^) is independent of Q, and therefore 
it suffices to consider the case O = 0, when the spheroidal wave functions reduce to the 
Legendre polynomials. Counting their zeros, we conclude that the function 

Y± has n zeros on (0, f ). (2.10) 



3 The Complex Riccati Equation 

Let Y\ and I2 be two real fundamental solutions of the Sturm-Liouville equation ()2.2|) for 
a general real and smooth potential V. Then their Wronskian 

w := Y l (u)Yi{u)-Y{(u)Y 2 {u) 

is a constant; we assume in what follows that w > 0. We combine the two real solutions 
to the complex function 

z = Y 1 + iY 2 , 
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and denote its polar decomposition by 

z = pe iv (3.1) 

with real functions p(u) > and f(u). By linearity, z is a complex solution of the Sturm- 
Liouville equation 

z" = V z . (3.2) 

Note that z has no zeros because at every u at least one of the fundamental solutions does 
not vanish. Thus the function y defined by 

y = - z (3.3) 

is smooth. Moreover, y satisfies the complex Riccati equation 

y' + y 2 = V. (3.4) 

The fact that the solutions of the complex Riccati equation are smooth will be helpful 
for getting estimates. Conversely, from a solution of the Riccati equation one obtains the 
corresponding solution of the Sturm-Liouville equation by integration, 

pv 

l °sz\'"u = y ■ (3.5) 

J u 

Using ()3.1|) in (|3.3|) gives separate equations for the amplitude and phase of z, 

p' = p Re y , if' = Im y , 

and integration gives 

V 



log p\ v u = Rey (3.6) 
J u 

< = flmv. (3-7) 



Furthermore, the Wronskian yields a simple algebraic relation between p and y. Namely, 
w can be expressed by w = Im (z z') = p 2 Im y and thus 



P- = • (3.8) 

im y 

Since p 2 and w are non-negative, we see that 

lmy(u) > for all u. (3-9) 

The boundary conditions for the Sturm-Liouville equation can easily be translated 
into conditions for y. To this end we write the solutions of the Sturm-Liouville equa- 
tion corresponding to the eigenvalues as = Im(e~' ia z^) with a suitable phase 
factor e~ ta . Then a Dirichlet condition can be written as <p = amodvr. A Neumann 
boundary condition gives 

= lm(e- ia yz) = Re (V i(Q+ 5 } y,z) 

7T 7T 

Rey cos((/j — a — — ) — Imy sm((p — a — — t 
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and thus 



7T 



a + — + arctan ( 



fRey 



\Im y 



According to ()3.7|) and ()3.9|) the function (p(u) is monotone increasing. Therefore, the 



number of zeros of Y, (|2.10j) tells us how often ip crossed the points mod it. This allows 
us to completely determine the "phase shifts" on the interval (0, ^), 



TT 



— + arctan ( 



Vim 



+ nil 



(n + l)ir 



(3.10) 
(3.11) 



(we use the usual convention that the arc tangent takes values in (—§,§))• Using (|3.7|) 
these boundary conditions can be expressed purely in terms of y and the integral of the 
imaginary part of y. 

For the gap estimates we need to control how y depends on A. To this end, we 
differentiate through the complex Riccati equation Q3.4|) and use that d\V = —1 according 
to (1231) and (J231). This gives the linear ODE 



y'x 



1 - 2yy A , 



where the A-derivative is denoted by a subscript. This equation can immediately be 
integrated using variation of constants. Applying (|3.5j) . we obtain 



z 2 y\\ 



Substituting the integration-by-parts formula 



we obtain the identity 



z 2 y\ 



„2y 



Z 

2y 



z 

2z7 



V 

+ 

u J u 



vy -y 2 z 2 

iy 2 



y-y 2 2 2 

iy 2 



In our estimates we will work both with (|3.12|) and (|3.13|) . 



(3.12) 



(3.13) 



4 Invariant Disk Estimates 

In this section we describe estimates for the complex Riccati equation (|3.4|) with initial 
conditions at u = 0, 

y = V-y 2 , y(0) = y (4.1) 

on the interval [0,u max ) with ii max G M + U {oo}. In what follows, the potential V £ 
C°°([0, «„,„)) can be any real and smooth function. The next lemma is the key to all the 
estimates in this section. 
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Lemma 4.1 Let a be a real function on [0,u max ) which is continuous and piecewise C . 
For a constant Tq > 1 we introduce the functions a, U and T by 



a{u) = exp I 2 J aj (4.2) 

U(u) = V-a 2 -a (4.3) 
T(u) = T exp Q TV [0)U) log\a 2 U\\ . 

Furthermore, we define the functions (3, R and m by 



(4.4) 



(3 = vM(r+±) (1.51 



a = m(T-l) (4.6) 



2 V T , 

m = a + ij3. (4.7) 

Suppose that U < on [0,u max ). If a solution y of the boundary value problem i4-l\ ) 
satisfies at u = the condition 

\y-m\ < R, (4.8) 
then this condition holds for all u £ [0, u max ) . 

Before coming to the proof, we briefly discuss the statement of this lemma. If a is a 
real solution of the Riccati equation, the function U as given by ()4.3|) vanishes identically, 
and thus = = R. In this case, the above lemma reduces to the trivial statement 
that y(0) = a implies that y = a on [0, n max ). It is more interesting to consider the case 
that a = Re y with y a complex solution of the Riccati equation. Then 

U = Re(V-a 2 -a) = Re (V - y 2 - y ) - (Imy) 2 = -(Imy) 2 < 0. 

Moreover, from (|3.6|) we can immediately compute a, 



a(u) = exp ( 2 / Rey 



o 



P 2 (u) 
P 2 (0) ' 



Substituting these relations into (|3.8|) yields 



a U = — itta ( Im y) " 



p 4 (o) v y/ p 4 (o) ■ 

Hence the function log \ p 2 U\ is a constant, and its total variation in 1)4. 4 J) vanishes. This 
means that T is a constant, and thus (3 and R are constant multiples of Imy. Our Lemma 
states that the circles of radius R{u) around the point m(u) = a{u) + i(3{u) are invariant 
under the flow of the Riccati equation. 

If no solution of the Riccati equation is known (and this will of course be the usual 
situation), one can take for a the real part of an approximate solution of the complex 
Riccati equation. In this case, the function log|/9 2 C/| will not be constant, but we can 
hope that its total variation is small. If this is the case, our lemma gives an "improved 
approximative solution" m together with a rigorous error estimate R. A good candidate for 
an approximate solution would be the usual wave function obtained by "gluing together" 
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suitable WKB wave functions and Airy functions as used in the semi-classical analysis of 
one-dimensional Schrodinger problems. We remark that the above lemma might even be 
useful for getting rigorous error estimates for numerical solutions of the Sturm-Liouville or 
Riccati equations. In this case, one would have to estimate the total variation of log \p 2 U\ 
from above, and this might be doable numerically if one has some control of the accuracy 
of the numerical calculation. 



ru 


\a 2 U\' 


Jo 


\a 2 U\ 



+ ee- u ) (4.9) 



Proof of Lemma \4-l\ For e > we set 

T e (u) = T exp Q 

and let R £ and m £ be the functions obtained from (|4.5jl — 1|4.7|) if one replaces T by T £ . 
Since T e (0) = T(0) and lim e ^oT e = T, it suffices to show that for all e > the following 
statement holds, 

\y - m £ \(0) < R £ (0) \y - m e \(u) < R £ {u) for all u G [0, u max ). 

In order to prove this statement, we will show that the assumption 

\y-m £ \(u) = R £ {u) (4.10) 

implies that 

\y-m £ \'{u) < R' £ (u) . (4.11) 

In what follows we will often omit the subscript e. 

Assume that Q4.1UJI holds and that U < 0. Then we can represent y as 

y = m + Re itp (4.12) 

with ip £ [0, 2ir). Furthermore, it follows immediately from (|4.5|) . (|4.6|) . and ()4.3jl that 

R 2 = U + 13 2 = V-a 2 + (3 2 -a. (4.13) 

Using the above relations together with (|3.4[) . we obtain 

- — \y — m\ 2 = (Re y — a) (Re y — a)' + (Im y — (3) (Im y — (3)' 



2 du 



Km 



(Rey-a) [V - (Re y) 2 + (Imy) 2 - a ] -(Imy-/5) [2 Re y Im y + [)'] 
(Re y — a)[V — (Re y) 2 - (Im y) 2 + 2/3 Im y - a] 
+ (Rey — a) 2(Imy — 0) Imy 

— (Im y — /3) [(3' + 2a Im y\ — (Imy — (3) 2(Rey — a) Imy 
(Rey - a) [V - (Rey - a) 2 - (Imy - [3) 2 - a 2 + (3 2 - a'] 
-(Imy - (3) [13' + 2a/J] - 2a ((Rey - a) 2 + (Imy - (3) 2 ) 

R cos 99 [V - R 2 + |m| 2 - a - 2a 2 ] - R sin <p [(3' + 2a(3\ - 2aR 2 

-2aR 2 -R{(3' + 2a(3) siny? < -2aR 2 + R\(3' + 2af3\ . 



Using that -j~|y — m| 2 = 2i?|y — m\' , we obtain the simple inequality 

\y-m\' < -2aR+ \0 + 2a(3\ . 
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Hence in order to prove (|4.11|) . it suffices to show that 

R' > -2aR + \(3' + 2a(3\ . 
Using ()4.2j) . we write the last inequality in the equivalent form 

(aR)' > \(aP)'\. 



(4.14) 



1 -(^/\aHT\T + ^aHAT- 1 



In order to prove this inequality, we first use (|4.5|) and (|4.6|) to write the functions a(3 and 
aR as 

a(3 
aR 

By definition of T e (JO), 



(4.15) 



\a 2 U\ 



\a 2 U\ 



+ ee 



It follows that 



(VWUIT- 1 )' 



-ee 



^(^laWlT- 1 ) 



(V\aW\T)' = se~ u (^\aHJ\T) 



if \a 2 U\' > 
if \a 2 U\' < . 



Hence when we differentiate through (|4.15f) and set e = 0, either the first or the second 
summand drop out in each equation , and we obtain {aR)' = \aj3\' . If e > 0, an inspection 
of the signs of the additional terms gives (J4.14|) . ■ 



The question arises how the function a in the above lemma is to be chosen. At this 
point, it is very helpful to regard as the one-dimensional Schrddinger equation for 

a quantum mechanical wave function Y, because this makes it possible to use ideas from 
semi-classical analysis. In order to explain our method, we first consider the WKB wave 
functions |H] 



i 



<f)(u) = |V|~s exp I ±i J y/\V\ 

which should be good approximations to fundamental solutions in the "semiclassical 
regime" V <C 0. The corresponding function y is 

&' V 

y(u) = — = ±i\/|V| . 

Lemma 14.11 should give a good estimate only if m is close to the exact solution y. This 
leads us to choose for the function a = Rem the corresponding expression in the WKB 
approximation, 

V 

a = in the "semiclassical regime" . 

4V 

This gives rise to the following estimate. 
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Theorem 4.2 Assume that the potential V is negative and monotone increasing on the 
interval [0, u max ), and that the following condition holds, 



K 



sup IF"! + TV U" V' 2 

h SUP < 1 



\v\ 3 



(4.16) 



where V max := supV < (and the supremum as well as the total variation are taken on 
the interval [0, u max ) ). Then the solution y of the boundary value problem with initial 
condition 

no) 



4V(0) 



satisfies on [0, u max ) the inequalities 



/ V 

y-*V\v\ + w 



Imy > 



< 20^\V\K 

W\ 
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Proof. We introduce on [0, u max ] the function a by 

V 

«(«) = 

Then from (|4~3"J) . 



o 



C7 



V" V 

+ 



12 



W AV 2 



Using the inequality Q4.16|) we get 



W 2 16|U| 3 



V 



2V < U < — 
~ 2 



In particular, U is negative. 

The inequalities (jOT}) and (jUHJ) allow us to estimate y/\U\ - ^/\V\, 



\U-V\ 



U-V 



V 



< 



< V\v\ 
'\u\ + V\v\ ~ 

Dividing by \f\V\ and \/|f7|, we obtain furthermore 

i fm 

i + K - y \u\ - 

Choosing Tq = 1 + K, we have the following estimates at u = 0, 



K 



\y — m \ 



\y\-p\ 
W\ 



{l + K) + 



1 



l + K 



l + K) + 



1 

l + K 




(4.17) 

(4.18) 
(4.19) 



(4.20) 
(4.21) 

(4.22) 



(4.23) 



(4.24) 
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Applying (|4.24l) we obtain 



\y — m\ < 



l + K) 



1 



l + K 



R. 



We conclude that the inequality ()4.8|) holds at u = 0. 

Hence we can apply Lemma 14.11 and obtain that (|4.8j) holds for all u G [0,u max ). 
Combining this with the inequalities 1)4.22(1 and (|4.24l) we obtain 

/I77T V 



< \ y -m\ + \p-y/\V\\ < R+\0-y/\V\\ 



R 



T + --2 




< 2y/\V\ I /' 
Jmy > (3-R 



1 



1 + if 



|c/| ^r- 

< 2y/\V\ {T -l + K) 




T ~ V2T ' 

It remains to estimate the function T, (|4.4j) . We first compute a and cr 2 ?7, 



(4.25) 
(4.26) 



a 



' Vo 
V(u) 



a 2 U 



V" 



5V 



rl 



| Vb| 4F 2 16|V| 3 ' 

where we set Vq = V(0). Applying (|4.16|) we immediately obtain the inequalities 



(4.27) 
(4.28) 



< 



a 2 U 



Vn 



< 2. 



The lower bound allows us to leave out the logarithm in the total variation in the definition 
of T; namely, 



TV [0tU) log\a 2 U\ 





(a 2 U)' 


< 2 


{a 2 U)> 


Jo 


a 2 U 


Jo 


\Vo\ 



2 TV 



a 2 U 



[o,uy 



Vn 



We substitute (|4.28|) into this equation and estimate the total variation of the individual 
terms using (|4.16j) as well as the monotonicity of V, 



V" 

TV [0,«)y2 ^ 



V 2 



+ 2 



u \ v"\ v Tvy + supiv"! 



{-Vf 



TV 



V 



[o,«) 



(\v\r 



< 



u 2\v"\v r3|y'i 3 supiy'i 

< — 1 - + 



(\v\)' 



v A 



V 2 



In the last term we can integrate by parts, 



u 3 \V'\ 3 
~V4~ 



v a ((-v)- 3 )' 



rl 



V 

1W 



f u 2V" V 



< sup — — 7 + 



u 3\V\ 3 
~V r ~ 



V' 2 sup | V" | 



y\Z 1 y2 
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Collecting all the terms and using (|4.16l) we conclude that 

TV [M log|a 2 £/| < 2K. 
We substitute this bound into (|4.4|) and use that T = 1 + K to obtain the bound 

T-l = (l + K)e K -l < \e K -l\+Ke K < 2e K . 
Using this bound in (|4.25f) and (|4.26[) concludes the proof. ■ 

The condition (|4. 16f) will clearly be violated when |V| becomes small. This is not 
astonishing because the WKB approximation fails near the zeros of the potential. In this 
"quantum regime", there is no canonical candidate for a, and therefore we simply take 

a = const in the "quantum regime" . 

We state the corresponding estimate in such a way that it can easily be "pasted together" 
with the result of Lemma 14.21 

Theorem 4.3 Assume that the potential V is negative and monotone (increasing or de- 
creasing) on [0, u max ), and that for some constant k > the following condition holds, 

VWo\u max < k (4.29) 

(with Vq = V(0)). Then any solution y of the boundary value problem \4-l\) which is 
bounded by 



1 2/0 1 < ci V | Vo| , Imy > 

ci 

for a suitable constant c\ > 1 satisfies on [0, u max ) the inequalities 

1 \V \ 



\y\ < c 2 Vploo , Imy > 



c 2 VWho 

where ||V||oo := sup[ 0)Mma:E ) |V| and C2 is a constant which depends only on k and c\. 

Proof. Let a be the constant function a = \J\Vq\. Then the function U = V — a 2 is 
clearly negative. A simple calculation shows that by choosing To = 2 ci(l + ci) 2 , we can 
arrange that |yo — m(0)| < R(0). Lemma l4~Tl yields that \y — m\ < R for all u G [0, it ma5c ). 
Since a is a constant, the function a is given by cr(u) = e 2au and thus 

\a 2 U\ = e 4au {a 2 - V) . 

As a consequence, 

\a 2 UY A V 

! < 4a 



\a 2 U\ ~ a 2 -V " 

If we integrate and use (|4.29|) . we obtain the following bound for T, 

T < T e 2au ( Vc^Vq Vh^V ] 



< rp Van ^VW 

a 
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Finally, we bound y by 

\y\ < \y — m \ + \ m \ < R + a + p 
= y/\U\T + a < (2T + l)a 



Imy > [5-R = > EL. 

These are the desired inequalities if we set C2 = 2T max + 1 = 8ci (1 + c\) 2 e 2K + 1. ■ 

It is obvious from (|2.3l 12.4(1 that the potential V has a singularity at u = 0. We now 
explain how Lemma 14.11 can be used for estimates near such a singular point. We will 
restrict attention to the case k = 0, but our method applies similarly to general k. In 
order to find a good candidate for the function a, we consider on the interval [0, u max ) 
the Sturm-Liouville equation with a potential which at u = ii max has the same singular 
behavior as 1)2.3(1 . 

1 

4 (t^max 

Setting u = u max — u, this differential equation has the two fundamental solutions = y/v 
and 4>2 = y/v logv, and therefore the function 



z"{u) = -— —^z. (4.30) 



z = \/v (1 — i log v ) 

is a complex solution. The corresponding solution of the complex Riccati equation is given 
by 

z' 1 i ( 1 logf 



^ z 2v + v (1 -ilogv) V 2u u (1 + log 2 v)J + u (1 + log 2 u) ' ^' 31 ^ 

Choosing a equal to the real part of this function gives rise to the following estimate. 

Lemma 4.4 Suppose that the potential V is on [0, u max ) of the form 

1 1 

with a function B which is monotone (decreasing or increasing) and satisfies the inequality 



V = --77. 7T2+ B (u) 



^LU + log'^HSlloo < I (4.32) 

(with ||-B||oo := su P[o,« moi ) Then any solution y of the boundary value problem fr4.1\ ) 

with initial conditions bounded by 

\yo\ < Cy/\Vo\, Imy > v ' 11 " 1 



C 

for any constant C > 1 satisfies on [0, u max ) the inequalities 

1,1 < -^r (4.33, 

Imy < 64 C 3 (1 + log 2 u max ) — — t (4.34) 

(u max - u) log [u„ 



1 



Imy > — 2 7 7 27 7- ( 4 - 35 ) 

64 C 6 (1 + log u max ) (u max - u) log (u v 



*- : rti t'l-r 
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Proof. We set v = u max — u and choose for a the real function 

1 log v 



a 



2v w(l + log 2 u) 



Using that a = Re y with y according to (|4.31|) and that y is a solution of the complex 
Riccati equation corresponding to the Sturm-Liouville equation ()4.3Uj) . we obtain 

U = V - a 2 - a' = Re (V - y 2 - y') - (Imy) 2 = B - — - 2 ^ 2 . (4.36) 

v z (1 + log v) z 

Using the assumption (|4.32|) together with the fact that the function v 2 (l + log 2 v ) is 
monotone increasing, we obtain that U is negative. 

At u = 0, the potentials V and U can easily be bounded from above and below, 

-- < -l-4v 2 \B\ < lAv 2 V = -l + 4v 2 \B\ < -- 
2 - i i - 2 

— — < Av 2 (l + log 2 v) 2 U = -1 + v 2 (l + log 2 v) 2 B < -- 
and in particular 



1 < (l + log 2 t;) y/\U\ < 

2 " IVUI ~ "' 



A simple calculation shows that by choosing To = 2C(1 + C) 2 (1 + log 2 n max ), we can 
arrange that |yo — m (0)| < R(0). Lemma l4~T1 yields that \y — m\ < R for all u E [0, u max ). 
Writing the function a in the form 



a = iL log (\A ^ + lo s 2u ) 

we can immediately compute a 2 U, 

a 2 = W 2 (l + log 2 t;) 2 
\a 2 U\ = 1 - v 2 (1 + log 2 u) 2 B . 

Using the bound (|4.32|) . we obtain 

TV [0iM) log|a 2 ?7| < 2TV [0 , n) |a 2 C/| < 4 n 2 ax (1 + log 2 u max ) 2 PIU < 2, (4.37) 

and thus T is bounded by T < T e 2 < 64C 3 (1 + log 2 n max ). Finally, we combine the above 
estimates with the inequalities 

\y\ < R+ \a\ +P , R-p < Imy < R + 0. u 

The estimate (|4.34jl is very useful because it shows that the pole of Imy at u = is 
integrable. 
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5 Convexity Estimates 

The estimates of the previous section gave us good control of the solutions of the boundary 
value problem (|4.1|) provided that the potential is negative. In this section we proceed 
with estimates in the case that V is positive, V > 0. Under this assumption, it is a simple 
observation that p 2 is convex, because 

(p 2 )" = {zzf = 2{V + \y\ 2 )p 2 > 0. (5.1) 

This fact will be essential for the estimates in this section. 
We begin with a lemma which bounds p from below. 

Lemma 5.1 Suppose that V is positive and monotone increasing on [0, u max ). Then every 
solution of the boundary value problem j4.1\) satisfies on [0, u max ) the inequality 

Imy 
P > Po — r 

w\ 

(with po = \z(0)\ and z any solution of the corresponding complex Sturm- Liouville equa- 
tion iTO) ). 

Proof. Differentiating the equation p' = p Re(y) and using the complex Riccati equa- 
tion (ITU) , we get 

p" = P (Rey) 2 + pRe(V-y 2 ) = (V + (Imy) 2 ) p , 
and using (|3.8jl we obtain the following differential equation for p, 

p" = Vp + ^. (5.2) 

r 

We let p{u) be the solution of the boundary value problem 

w 2 

p" = ^ with p(0)=p , p'(0) =p'(0) = Po Rey . (5.3) 
-pi - 

The function p is a solution of (|5.2|) in the case V = 0. Therefore, p can be written 
explicitly in the form p = \z\ with z a solution of the complex Sturm-Liouville equation 
without potential with Wronskian equal to w, i.e. 

z" = and Im(|z') = w = p^Intyo . 

A short calculation shows that p has the simple form 

p(u) = po \l + y u\. 

This function is defined even for all u G R. It is convex, and its minimum is computed to 
be 

mm p(u) = po -, — r • (5-4) 
ueR- \y \ 

We introduce the set 7 C K 2 by 

I = {(p{x),[f/(x),oo)) with x G R} . 
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Let us show that I is an invariant region for p in phase space, i.e. that for all u G [0, u max ), 

(p(u),p'(u))el. (5.5) 

Once this is shown, the lemma follows immediately from ()5.4|) . Due to our initial condition, 
(|5.5|) is clearly satisfied at u = 0. Thus assume that [0, v ] with < v < u max is the maximal 
interval where (|5.5j) holds. Then the point (p(v),p'(v)) lies on the boundary of /. Using 
that p is convex and thus p' is monotone increasing, one finds that either 

p(v) = minp and p'(v) > 

or else there is i 6 R such that 

p(v) = p{x) and p(v) = p'(x) . 

In the first case, it is obvious that the gradient of (p(v),p'(v)) is pointed towards the 
interior of /. In the second case, the estimate 

11/ \ T 7" / \ W W W II I \ 

p (y) = Vp(v) + —-^ > —-3 = -—3 = p (x) 

Plv) p{v) P\ x ) ~ 

yields that 

p (v) = p'(x) and p"(v) > p"(%) ■ 

Hence the gradient of (p(u),p'(u)) again points towards the interior of /. We conclude 
that ()5.5j) holds also in an interval [v, v + e) with e > 0, a contradiction. ■ 

This lemma has the following immediate consequence. Due to the convexity of p, 
supp < p + p{u) = p(u) ( 1 + -7^- J < p(u) ( 1 + //n 



[0,u) V P( U )J \ lm Vo 

and hence 

sup p < p(u) —^2- . (5.6) 

io,u) Im yo 

In regions where the potential V is large, we expect that p should increase exponen- 
tially. The next lemma quantifies this exponential increase of p by showing that in the 
"semiclassical regime" V 3> 0, the integral over p 2 is much smaller than the supremum 
ofp 2 . 

Lemma 5.2 Suppose that V is positive and monotone increasing on [0, u max ). Then every 
solution of the boundary value problem \4-l\ l satisfies on [0, u max ) the inequality 

p 2 < L sup p 2 

[0,w) 

with L given by 

( 3 V'\ V 
L = sup^ + ^j+TV^. (5.7) 
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Proof. We substitute the differential equation for p 2 , (|5.1|) . into the integral, 



1 



1 



Integrating by parts gives 



(p 



2\l 



l{V + \y\ 2 ) 



2 7 V + \v\ 



(p 



2\n 



1 



v + \y\ 2 



(p 2 y ■ 



Using the estimates 

(p 2 y 



V+\y\ 2 
1 



V + 1 2/| 2 



< 



< 



2p 2 \y\ < 2p 2 |y| 



^ + |y| 2 ' : 2VV\y\ Vv 

V' + 2\y\ \V-y 2 \ V 2\y\ 



(V + \y\ 2 ) 2 



< 



+ 



v l 

< ^ + 



V 2 V + \y\ 2 ~ V 2 VV 



we obtain 



p 2 i p ( V 1 



+ 



+ 



2\l\ 



When integrating by parts once again we must be careful because the function (p 2 )' may 
change signs. However, since p 2 is convex, it changes signs at most once, and therefore we 
get positive boundary terms at most twice, 



u (V l 

\v 2 ~ + Vv 



d_ 

du 1 



[0,u) 



v 2 Vv 



+ 



U 2 + Vv 



Finally, we can estimate the last integral by 



Y_ l 

v 2 + Vv 



p 2 < 



2 TV (V , 1 

suppTV [M ^ + w 

2 L v , i 

supp TV [0u )-2 + sup^= 
[0,u) V ^ [0,u) W 



where in the last step we used the monotonicity of V. 



6 Elementary Properties of the Potential 

In this section we shall analyze the potential V (|2.3M2.4j) for large A and f2. More precisely, 
we consider the range 

n > n and A > 2Afl (6.1) 

for parameters fio an d A, which we can choose as large as we want. Then the potential 
looks qualitatively as in Figure ^ In the case k ^ 0, V has a unique minimum uq given 

by 

sin 2 Uo = 1^2-1, (6.2) 
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and the potential is negative at the minimum, 

(6.3) 




(where in the last step we possibly increased A). V is strictly decreasing on the interval 
[0, uq] and tends to infinity as u \ 0. Thus there is a unique u_ G [0, uq] with V(u-) = 0. 
On the interval [no, § ], V is strictly increasing. Thus there is at most one u + G (no, ^] 
with V(u+) = 0. If no such u+ exists, we set u + = ^. For a given parameter k > 
(which will be specified later) we set An = k/V 'Ail. It is easily verified that by choosing A 
sufficiently large we can arrange that V(uq ± Au) < 0. As a consequence, \u± — uo\ > Au 
and thus, using (|6.3|) . 

|V(uo)| |u± - u | 2 > k 2 . 

Using monotonicity, we can thus uniquely introduce points u+ G (uo,n+) and u_ G 
(n_, no) by the condition that 

\V(ui)\ \u± - u|| 2 = k 2 . (6.4) 

Finally, we introduce the point u 1 G (n+,§] by the condition Viu 1 ) = Oa. If no such 
point exists, we set u 1 = ^. 

In the case k = 0, V is monotone increasing on the whole interval (0, We set 

j 1 k 1 

n = — = — , i*o = —p^ and u\ = — = . (6.5) 

8VA log 2 A V\ 

The points n + , u+ G (ni,n + ) and n 7 ^ are introduced as in the case k 7^ 0. 

We consider on (0, |] the solution y of the complex Riccati equation (|3.4|) with initial 
condition 

y(«o) = ^vT^oI-^y- (6-6) 

The next lemmas make the following statements precise: The intervals S (as introduced 
in Figure ^) are "semiclassical" in the sense that Theorem 14.21 applies. On C we can 
use the convexity estimate of Lemma 15.21 The interval P near the pole can be treated 
by Theorem 14.41 Finally, the "intermediate regions" I± and J can be controlled with 
Theorem 14.31 and Lemma 15. II 

Lemma 6.1 For every 5 > and k G Z there are parameters k, A, Qq > such that for 
all SI, A in the range \b\l\) . the quantity K as defined by \4-lb\) is on the interval S bounded 
by 

K < 5. 



Proof. The third derivative of V can be written in the form 

cos u 
sin 5 u 



V'" = —t~k — (polynomial in sin 2 u of degree 3) . 



Hence V" has on any interval [u, v] C [0, 5] at most 4 zeros. Thus, after splitting up 
[u, v] into at most four subintervals, V" has on each subinterval a fixed sign. On any such 
subinterval \u, v\ we can apply the estimate 



[ V \V'"\du < \V"(u)\ + \V"(v)\. 

J u 
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This makes it possible to control the total variation of V" in (|4.16j) by 8sup|U"|. We 
conclude that it suffices to show that on the interval / the following two inequalities hold, 

< S (6.7) 



\V\ 3 
W"\ 



< 5. (6.8) 



V 2 

We treat three cases separately. 
First case: k = and u + > 

On the interval [||-, the potential V is concave; more precisely, 

-n 2 < V" < -—. (6.9) 
4 

Integration yields for all r £ (j^, u+] the following bounds for V and V, 
£l 2 fir \ rrl/ , . 2 



^r(^-r) < V'(r) < 2 (~-r) (6.10) 

n 2 



Q 2 /7T \2 



Since V(u+) is either zero or negative, it follows that 



(u+ — t) {-k — u + — t) 



\V{t)\ > ^-( u+ -r)(7T-u + -r). (6.11) 
Combining the inequalities (|6.91 16.101 16.11jl . we obtain for all r £ (t§" » the estimates 
F'(r) 2 . ft 4 (f-r N2 



|F(r)| 3 " ^ 6 («+-t) 3 (7t-« + -t) 3 " ^2( n+ _ u 5 )4 ^' ' 

\V"(r)\ & 64 

F 2 (r) " tt 4 (u+ - r) 2 (vr - u + - rf ~ Q 2 ( u+ - u s + )* ' [ ' 

In order to estimate the factor (u+ — u+) from below, we use 1)6.11)1 in the defining equation 
for u'l, (fOj) . 

^(u + -u s + r < \V(u s + )\(u + -u s + ) 2 = k 2 . 

Using this inequality in ()6.121 16.131) and choosing k sufficiently large, we obtain 1)6.71 16.8)1 
for all r £ (ff,^]. 

On the interval [uo,Ui], a short calculation using (|2.31 16.51 16.1)) shows that 

A A 3 2A 2 

\V\>~, \V'\ 2 <-, \V"\ < — on[u ,ux], (6.14) 

A Kj K/ 

again proving ()6.71 16.8)1 . 

i +■ an ra 1 I ij-i 

16 



On the remaining interval (ui, 15), we know from the monotonicity of the potential 



and (KTTT1) that 

Q 2 ( 

XV(r\\ >_(«,_ , , , . 

16 7 V 16 7 _ 2048 



|V(r)| > — «+-— 7T-U+-— >t^0 2 
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Furthermore, a short calculation using (|2.3|) shows that on (u\, ^) 

\V'\ + \V"\ < 40 2 on [m,f] . (6.15) 

We conclude that, choosing Qq sufficiently large, we can again arrange that (|6,71 I6.8|) 
holds. 

Second case: k = and u + < 

On the interval [«o,ui], we can again use the estimate (|6.14|) . Conversely, on the inter- 
val [tii, xh a snor t calculation shows that V" can be bounded in terms of higher powers 
of the first derivatives; more precisely, 

\V"\ < 10 |V'|! . 

This inequality allows us to deduce (|6.8|) from ()6.7|) . Hence it remains to prove the in- 
equality (|6.7|) on the interval 

On the interval the potential V is either convex or else at least the second 

4 

derivative of V is large compared to \V |s. More precisely, by choosing Qo sufficiently 
large, we can arrange that 

V" > V'* . (6.16) 

We shall derive an upper bound for Aw := ii+ — u+; for ease in notation the subscript '+' 
will be omitted. We rewrite (|6,16|) as 



du\ ) 3 y'% 

We integrate from u s to r G [n 5 ,u] to obtain 



d / , i\ 1 V" K'i 

< 



_2 \ -3 

K 3 



V(t) > [V'(u s )^+^(r-u s )j . 

Integrating r over the interval [it s ,u], we obtain for AV := V(u) — V(u s ) the estimate 

- / \ — - 

AV > — V\u s )l [ 1 - , 1 \ with a := — ^V) 1 An . 
~ 2 v ; V (1 + a) 2 J 3 V 7 

The inequality 

i-— ^> Q 



(1 + a) 2 ~ 1 + a 
gives 

AV > V ' {US) AU 



2 1 + a 

By definition of u s , (|6.4[) . we know that Ay -(An) 2 = k 2 . Hence, multiplying the last 
inequality by (Au) 2 , we obtain 

2 UV) (An) 3 
2 1 + a 

Using the definition of a gives the inequality 

(Au) 3 - ^ (J V'(u s )-^ 2 Au-2^J V'(u s )-^y < 0. 



23 



Since the polynomial x 3 — ^ — 2 is positive for x > 2, we conclude that 

Au < 2 k§ F'(-u 5 )~3 . 
Using again the relation AV(Au) 2 = k 2 , we get an upper bound for AV, 

2 

AV > ^V'{u s )i . (6.17) 

This proves the inequality (|6.7j) at u = u s . 

Next we shall show that Q6.7|) holds on the whole interval (u\,u s ]. To this end, we 
introduce on this interval the function / by 

4 3 

/ = V' 2 + V 3 . 

K 

We saw above that f(u s ) < 0; our goal is to show that / < on (ui,u ]. Let (v,u s ] 
with u\ < v < u s be the maximal interval on which / is negative. We apply ()6.16[) to 
obtain 



, ( 2 ,4 3 • 4 3 a\ 2,7 

> V I -2k-zV'z + — — V 2 \ > Wk-zV'3 > 0, 

where in the last line we used that f(v) < 0. The last inequality contradicts the maximal- 
ity of the interval (v,u s ] unless v = u\. This concludes the proof in the second case. 

Third case: k ^ 0. 

On the interval (u_ , uo] , the proof of the second case goes through without changes if 
we consider the integral backwards and set u\ = uq — l/(4vT2). On the remaining 
interval [iio,ti+], we can use the proof of the first case and the second case after set- 
ting m = u + l/(4\/n). ■ 



Lemma 6.2 For sufficiently large A and Qq, the parameter L (as defined by j5.7\) ) is on 
the interval C for allQ,, A in the range \6.1\) bounded by 




Proof. Similar to (|6.9M6.1(?|) . one easily sees that V' and V" satisfy on C the bound \V'\ + 
\V"\ < AVi 2 . On the other hand, it is clear from the definition of u 1 that |V| > fl? on C. 
This immediately gives the lemma. ■ 



Lemma 6.3 For sufficiently large A and Oo, the potential 

B{u) := V{u) + ^ 
satisfies on the interval P the inequality 



(u J ) 2 (l+log 2 u J ) 2 \\B\\ oc < I. (6.18) 
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Proof. A short calculation shows that B is bounded from above by 

\B(u)\ < n 2 u 2 + 2X. 
From the definition of u J , l|6.5jl . it is clear that for large A, | logn J | < log A, and thus 

(i* J ) 2 (1 + log 2 ^) 2 HSIloo < ) (1 + log 2 A) 2 ||g||oo 

64 A log A 

i „ „ Em wo 2 \ i 

< — r \\ B Woo < — r < 

~~ 32 A ~ 32 A V A ~ 

where in the last step we again used (|6,1|) . 



Lemma 6.4 For every 5 > and k £ "L i/iere are parameters k, A, Oo > snc/i i/iai /or 
aZZ O, A in iae range \6.1)) . 

|/±| < 5- 

Proof. We choose Oo so large that no < f • Then clearly |7_| < 5. Furthermore, it 
is readily verified that the potential is increasing on the interval K := [f , ^ — |] at the 
following rate, 

V(v) - V(u) > c (v - n) 2 Q 2 for all n, v £ K, v > u , 
where c is independent of A and 0. This implies that 

because otherwise (|6.4|) would be violated. Furthermore, the condition V(u J ) < O2 implies 
that 

|[n + ,u I ] n K\ < c~i f}~3 . 
We conclude that by increasing Oq, we can arrange that \I+(lK | < | and thus < 5. ■ 




7 Spectral Estimates for the Selfadjoint Problem 

In this section we shall prove Lemma f2.ll We begin by reducing the problem to an estimate 
for y x . 

Lemma 7.1 Assume that for any given k £ Z and e > 0, there are constants A, VLq > 
snc/t that 

/ 2 Imy A < £ (7.1) 
Jo 

for all £1 and A in the range \6.1\) . Then Lemma \2. 1\ holds. 
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Proof. According to the asymptotics (|2.9|) it suffices to consider A in the range A > AQ 
for sufficiently large A. Let us consider the quotient Re(y)/Im(y) in (|3,1U|) , According 
to Theorem 14.21 and Theorem 14.31 there is a constant c > such that for all f2 > J7q 
and A > AfJ, 

Rey 



Imy 

2 



> —c on [uq,u + ]. (7.2) 



On the interval [u+, f], p 2 is convex, and using the identity 



Im y w 2w 

one sees that Re(y)/Im(y) is monotone increasing. We conclude that the inequality (|7.2|) 
also holds at u = f , and thus 



7T 



+ + /Rey(f)\ 

< — arctan c < arctan < — . 

2 VImy(f)y 2 

Using the last bounds in (|3.1()j) one sees that for two neighboring eigenvalues, the 
phases must differ at least by 5 := ^ — arctan c, 

7T 

2 > 8. (7.3) 



From ()3.11j) one sees that this inequality is also true for the states of odd parity (with 5 = 
7r). Applying ()3.7|) and the mean value theorem, we conclude that there is A G [A^, A^ +1 
such that 



(A± +1 -A±) / Imy A > 5. 
Jo 

Hence the upper bound H7.1|) gives the desired gap estimate. 



2 



We establish (|7.1|) by deriving separate estimates in the regions S, I±, C and near the 
pole. 

Lemma 7.2 For any given k G Z and e > 0, there are constants A, 0,q > such that on 
the interval S, 

\y\\ < £■ 

Proof. Differentiating the initial condition (jfi.flj) gives 



+ 



V'(u ) 



This can be estimated using Lemma Id. 11 



\yx(u Q )\ < 



m^o)i 



i + 



< 



For given u £ S, we compute y\{u) via (|3.13|) . This gives rise to the estimate 



\V\( U )\ < 



P 2 (up) 
p 2 (u) 



y\(uo) + 



1 



^ 2Imy(u) I 



|V(«o)| 2|y(«)| 



2y(^o) 

+ Imy (it) 



1 1 

+ 



2|y(n)| (?{u) 

u \v-y 2 \ 



u \V -y 2 \ 
uo 2\yf 
1 



2 |y| 2 Imi/(n) 



(7.4) 



(7.5) 
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\y\ < 22 J\V 



where in the last step we used (|7.41 16.41 13. 8|) . According to Lemma 16.11 we can apply 
Theorem 14.21 We thus obtain for any K G (0, 1) the estimates 

(7.6) 
(7.7) 



(7.8) 



y-iy/W\ < \A^\ (20K + \^K) < 2lVKy/\V\ 

\V-y 2 \ = \{y - iVW\){y + iVW\) 

< 21^K^/\V\23 ^/\V\ < 500 |V| VK , 

and, after choosing K < 1/42, the inequality Q7.7JI shows that 

Imw > ^ ■ . 
y - 2 

Using the inequalities (|7.6l I7.9|) in (j7.5|) gives 

1 



I^AWl < — rrb — h 



+ 22^/\V{u)\ 



i^mi VWW\ 

and, since V is monotonous, it follows using (J7.91 17.8|) that 



u \V-y 2 



\y\( u )\ < 



< 



45 f u \V-y 2 
+ 22 - — 

VWM Juo \y\ 2 

45 J— 7T 

; + 88 • 500 - 



(7.9) 



The last expression can be made arbitrarily small according to Lemma 16. II 



Lemma 7.3 For any given k £ Z and e > 0, there are constants A, £Iq > smc/i that on 
the intervals I±, 

\y\\ < £■ 

Proof. We only consider the interval I+; the proof for I is analogous. For any v £ 
(u+,u+], we compute y\ via (|3.12|) with u = u+, 

»w = — ^ i^jX^ • 

According to the definition of u+, (|6.4jl . we can apply Theorem l4.31 This gives the estimate 

\yx(v)\ < \y\(u S + )\ Imy(V ] +lmy(v) f ^— < c| \y x (u% )| + c| (v - uf) . 
Imy(u%) J u s^ hay 

This can be made arbitrarily small according to Lemma 17.21 and Lemma 16.41 

If v € (u + , u 1 ], the change of y\ on the interval (u+,v) can be estimated similarly using 
Lemma 15.11 ■ 
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Lemma 7.4 For any given fcsZ and e > 0, there are constants A,Qq > such that on 
the interval C , 

\y\\ < e • 



Proof. We again compute y\ via ([3.12ft . This gives for any v G C the estimate 

P 2 (V) P 2 {V) J U I 

The first summand can be made arbitrarily small according to Lemma l7.3l and Lemma [5. II 
whereas the second summand can be handled with Lemma 15.21 and Lemma 16.21 ■ 



Lemma 7.5 For any given k G 7L and e > 0, there are constants A, f2 > such that for 
all A and £1 in the range \6.1)) . 

/ \v\\ < e - 

Proof. A short calculation using (|2.31 I6.5JI shows that on the interval J, 

^2 ^ l^(«o)l ^ 8 A log 4 A < \V(u J )\ < 64 A log 4 A, | J\ < . 

In particular, |l^(«o)| |J| 2 < 1, and so we can apply Theorem 14.31 to obtain on J the 
estimates 

\y\ < 16c 2 V\ log 2 A , Imj/ > — — ^- . (7.10) 

16c2/^ log A 



These estimates allow us to bound y\ on J again using (|3.12|) . Namely, for all v E J, 

Imy(v) / ,Uo 1 



Imi/(f) Z" 1 

< \v\{u )\- — -+lm.y{v) / 

lmy(u ) J v 



Im y 



Estimating the factor \y\(uo)\ by (|7.4|) . we obtain 



|y A («)| < on J (7.11) 

with C3 = 256c2/t 3 . By increasing A this can be made arbitrarily small. 

On the interval P, we apply Lemma f4.4l with C = 64c2/-c 2 log 4 A. This gives the estimate 

— ; — s - ; 0— < i m < 9 — ( ' -12) 

c 4 log p A v log v v log u 

with p = 14 and a constant C4 which is independent of A and 0. We again estimate y\ 
using (EHJ), 



Iva(«)| < \yx(u J )\ 



Imy(v) 



Im^u- 7 ) 

Estimating y(it J ) and y\(u J ) by (|7.10[ l7~TT|) and using (|7.12|) we get for all v 6 P, 

lo g p+6 A 1 log 2p A Z 1 "' 7 

\ i 2 + Cs 1 2 

A v log w w log t> 



< c 5 1 2 + c 5 ,_ 2 / T log 2 Tdr 
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for a suitable constant C5. This expression is not bounded as v \ 0. But the pole is 
integrable, and the calculation 



,J dv 1 







v log 2 v \ogw J 



l '' dv f uJ f uJ f T dv 

— — 2— / T !og 2 rdr = / dr r log 2 r / — — 5— 
f log v J v Jo Jo v log V 

= - t log rdr = - (u J f (1 - 2 log it J ) 
Jo 4 

shows that, by increasing A, we can make the resulting integrals over P arbitrarily small. 



This completes the proof of Theorem 11.31 

8 Slightly Non-Selfadjoint Perturbations 

It remains to prove Theorem ll.il In preparation, we split up the spheroidal wave operator 
as 

A = Aq + W 

with 

A = --r^— a sin 2 ??-^— + —L-(R e nsm 2 i} + k) 2 
d cos v d cos v sin v 

W = 2iImfi(ReOsin 2 i? + /c) - (Im S!) 2 sin 2 . 

The symmetric operator _4o clearly satisfies the hypothesis of Theorems I1.2I and II .31 
whereas the complex potential W is uniformly bounded according to assumption (|1.5|) . 

\W\ < 2(A; + l)c + c 2 =: |. (8.1) 

Our method is to treat If as a slightly non-selfadjoint perturbation as introduced by 
Kato 8, V.4.5]; see in particular Theorem 4.15a]. Unfortunately, the latter theorem 
is not quite strong enough for our purpose. For clarity, we here repeat the basic ideas of 
Kato and give a detailed proof of our main theorem. 

Proof of Theorem li.il Throughout the proof, we restrict all operators either to H+ 
or 7i-. Applying Theorems II .21 and II .31 to the operator Ao and 7 = 8p, we obtain for the 
eigenvalues < Ai < A2 < • • • of Ao the estimates 

A n +i — A n > 7 for all n > N and flgK. 

For all A a(Ao), the resolvent := (A — -4o) -1 exists and satisfies the bound 

llR ° A ~ dist(A,l(A)) ' (8 ' 2) 
Since the spectrum of Ao is real, we have in particular 

WR^W < ' 



ImA| 
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Around each A n , we draw a circle of radius p. The first N circles may intersect, and 
we take the outermost lines to define the contour Co, 

C = d(B p {\ x ) U . . . U B p (X N ) . 

All the following circles do not intersect and give rise to the contours 

C k = dB p (\ N+k ) , k>l 

(see Figure HJ. Since the distance of these contours to the spectral points of ^4o is at 
least p, we have for A on any of these contours, 



\R\W\\ < \\R° X \\ \\W\\ < 



Hence the operator 1 + R^W is invertible with a Neumann series. We conclude that the 
resolvent R x := (A - A)- 1 = ((A - A ) (1 - B^W))' 1 = (1 - R° X W)- 1 R° X exists and 

II^aII < 2||i?°||. (8.3) 

This allows us to introduce the operators Q k as the following contour integrals 

Qk = 7T— </ R\d\. 



27ri Jc k 

The Cauchy integral formula together with the resolvent identity 

R\ R\> = — ^ _ y (R\ - R\>) 

immediately yield that the operators Qk are projectors onto invariant subspaces of A, 
and that they are holomorphic in £1. Furthermore, they are uniformly bounded because 
according to (|8.31 18.2|) and the definition of the contours, 

\\Qk\\ < 7^ / 2||i?°|| < 2N. 

We introduce the operators Pk as the finite sums 

K 

Pk = Y^Qk- 

fc=0 
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X 



N+K 



Ai 



R 



Figure 3: The contours Dk- 



For the unperturbed operator Ao, we introduce similarly the projectors and P K . Let 
us derive estimates for the difference Pk — Pk- We first write it as the contour integral 

Pk~P°k = (R\~ R°x) dX , (8.4) 

2m Jd k 

where Dk is a rectangle with side lengths Xn+k + Xn+k+i and 2R centered at the origin 
(see Figure |2J). Since dist(Dfc, a(Ao)) > p, the inequality ()8.3|) again holds. Using the 
resolvent identity 

R x - R° x = R x W R° x 

together with (|8.1( I8.2j) , we obtain for any A on the contour Dk for sufficiently large R 
the estimate 

\\R\ - Rx\\ ^ 



- 7 2 + 4 (ImA) 2 ' 

This inequality allows us to take in (|8.4|) the limit R — > oo to obtain the estimate 

\\Pk-P° k \\ < ^ = \- (8.5) 

This estimate can be improved if the operator Pk — P K is restricted to the range of P°, 
L < K. Namely, applying the bound 

\\R° x Pl\\ < max IA-AJ" 1 

to the resolvent identity gives for any A on the contour Dk for sufficiently large R the 
bound 

11(^-^)^11 < | ^ ? . T 8 " T r. 

|7 + 2z ImA | |A - Aat + l| 

Substituting this estimate into the contour integral, taking the limit R — > oo and estimat- 
ing the resulting integral as follows, 

dx dx f°° dx 2 



Va 2 + x 2 Vb 2 + x 2 Jo ab x 2 yfab 



we conclude that 



\\(Pk ~ P° K ) P°l\\ < , n 8P . ^ , i < K. (8.6) 
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The inequality ((8.5|) allows us to determine the rank of the operators Pr. Namely, for 
every * in the range of P K , 

\\p K *w > \\p^\\ - \\(p K -p° K )n > i n*n . 

In particular, * is not in the kernel of Pk- This shows that the rank of Pk is greater or 
equal to the rank of P K . Interchanging the roles of Pjc and P K , we see that Pk and P K 
have the same rank. Since P K is the projector on the eigenspaces of *4o corresponding to 
the eigenvalues Ai, . . . , \n+k, the dimension of its range is N + k. We conclude that Qo 
is a projector on an iV-dimensional invariant subspace of A and Q\, Q%, . . . are projectors 
on 1-dimensional eigenspaces. 

The inequalities (|8.5U8.6|) imply completeness: Let * G Ti and e > 0. Since the spectral 
projectors of the unperturbed problem converge strongly (i.e. s — lirn^oo PjJ = 1), there 
is L G N such that ||* - < s for all M > L. According to (|8Jt|) . \\P k - 1|| < 

\\ p K ~ P K \\ + \\P K \\ + ||1|| < 3. Hence for all K > L, 

\\{p K - 1) *n < \\{p K - 1) (* - pin\ + \\{pk - p° K ) pin 

< 2e + ||(P/c-P&)Pg*||, 

and the estimate (|8.6|) shows that the last term can be made arbitrarily small by choos- 
ing K sufficiently large. 

It remains to prove the inequalities ()1.6|) and (jl.7|) . Combining the gap estimates of 
Theorem 11.31 with the fact that G B p (Xj^ + i t ), one immediately obtains Q1.7JI . The in- 
equality (| 1.7(1 follows similarly from the asymptotics ()2.9jl . ■ 
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